Cosmological parameter estimation with large scale structure and supernovae data 
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Most cosmological parameter estimations are based on the same set of observations and are there- 
fore not independent. Here, we test the consistency of parameter estimations using a combination 
of large-scale structure and supernovae data, without cosmic microwave background (CMB) data. 
We combine observations from the IRAS 1.2 Jy and Las Campanas redshift surveys, galaxy peculiar 
velocities and measurements of type la supernovae to obtain h = 0.571q'i4i fi m = 0.28 ± 0.05 and 
as = 0.87^0 05 in agreement with the constraints from observations of the CMB anisotropics by the 
WMAP satellite. We also compare results from different subsets of data in order to investigate the 
effect of priors and residual errors in the data. We find that some parameters are consistently well 
constrained whereas others are consistently ill-determined, or even yield poorly consistent results, 
thereby illustrating the importance of priors and data contributions. 



I. INTRODUCTION 

The measurements of the cosmic microwave background (CMB) from the WMAP satellite pj and large-scale galaxy 
clustering from the 2dF 0] and SLOAN surveys have provided tight constraints on the value of cosmological 
parameters The combination of observations of different nature has the power to tighten constraints on the 

parameters of cosmological models (see e.g. 0], H^, HI). A shortcoming of this approach, however, is 

that most parameter estimations use the same data. They are not independent and their agreement comes as no 
surprise. Also, if one set of data suffers from uncontrolled systematic errors, those errors might affect the results and 
go unnoticed. It is therefore useful to check those results with an independent set of observations. 

Some cosmological parameters such as the Hubble parameter h are ill-determined from one type of observations 
alone and a joint analysis is necessary in order to determine its value. If such an analysis is not available some prior 
value is often assumed, based on some dedicated observation. For example, the HST key project constraint on the 
value of the Hubble constant Hq = 0.72 ± O.O8/1 km s _1 Mpc _1 often helps to tighten the determination of other 
parameter values [2|. The robustness of such results might depend strongly on the prior. 

In this brief communication we use cosmological data to estimate cosmological parameters and to provide an 
independent test of more recent analyses. We combine data from the IRAS 1.2 Jy [2(J and Las Campanas p| redshift 
surveys, galaxy peculiar velocity measurements and measurements of type la supernovae 0, |J. We do not use 
any cosmic microwave background data. We examine our results according to subsets of data. We focus on three 
parameters: the matter density £l m , the Hubble parameter h and the amplitude of mass density fluctuations ag. 

Our paper is organized as follows. In the sections HTl and ITTT1 we describe the data and method used in the analysis. 
The results are presented in section ITv! and conclusions are presented in section Ivl 

II. LARGE-SCALE STRUCTURE AND SUPERNOVAE DATA 

The IRAS 1.2 Jy data set is an all-sky (87.6%) infra-red survey of 5313 galaxies. The density of galaxies in three 
dimensional space is sensitive to the mass density fluctuation parameter as and the inferred matter power spectrum 
depends on the shape parameter T. The likelihood function was calculated for flat universes within linear theory by 
20] and assumes that light traces mass via a linear bias factor (8p/p) = &iras (?>p/p) m - The likelihood is a function 
of the 3 parameters erg, /3iras and r. 

The Las Campanas Redshift Survey (LCRS) probes the same parameters as the IRAS survey but both the data 
and the likelihood function calculation differ in several ways. IRAS detected galaxies in the infra-red whereas the Las 
Campanas experiment measured ~ 26,000 galaxies at optical wavelengths [13j. A different morphological selection 
took place in each of the experiments such that their bias parameters are different. The LCRS probes a deeper 
range of redshifts in finger-shaped volumes consisting of a northern and a southern sample. A Karhunen-Loeve (KL) 
eigenmode basis was constructed to correlate the measured redshifts as described in Q (see also ^J)- The data 
were compared to matter power spectra expressed in terms of the KL basis from numerically simulated open-CDM 
structure formation scenarios with different values of ag, Plcrs and T. 
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Unlike redshift surveys which measure the distribution of visible matter, peculiar velocity measurements are a 
dynamical probe of the matter power spectrum. Within the assumption that galaxies trace the large-scale gravitation 
induced velocity field it is expected that velocities obey linear theory better than density due to the non-local character 
of the gravitational field. This likelihood function was calculated in using the Mark III catalogue of more than 
3000 galaxy velocities with no morphological selection (see |2l| and references therein) . 

The SNe la likelihood function used is calculated in and by maximising the likelihood function using 42 
high redshift supernovae measured by the Supernova Cosmology Project ^lj. Note that SNe la measurements 
probe cosmological parameters through the luminosity distance. Unlike redshift surveys, SNe la are independent 
of the structure formation scenario so the SNe la and LSS combination of data is a powerful tool for constraining 
cosmological models. 

III. METHOD 

In order to perform a joint analysis, one must ensure a basis set of independent parameters from which the other 
parameters can be calculated so the parameter space depends on the combination of data sets. 

The IRAS likelihood function was computed under the assumption of a flat universe. Hence only flat models are 
considered when IRAS data are combined with, e.g. supernovae data even if the latter likelihood function allows 
more freedom in the parameter space. Similarly, the Las Campanas likelihood function was calculated by comparing 
realisations of the survey from cold dark matter structure formation simulations. Therefore no hot or warm dark 
matter models can be constrained with the Las Campanas likelihood function. 

The assumptions of each data set are propagated throughout the analysis because they change the way the parameter 
space can be defined and how the remaining parameters are calculated from the independent ones. This conservative 
approach in which only models common to all used data sets can be constrained, has the drawback that the results 
cannot always be compared at face value. The limitation to flat universes corresponds to a slice of the larger parameter 
space allowing for open and closed universes. Results from the larger parameter space, marginalised over curvature, 
cannot be considered statistically equivalent to results for flat universes only. The priors on the considered models 
for each combination are summarised in table [I] 

If the likelihood functions are not evaluated beyond some value the nearest border value of the likelihood function 
is used. Although the IRAS likelihood function is incomplete in the eg dimension, this choice of extrapolation has no 
consequences when using the full combination of data sets. The analysis is carried out with Bayesys, a commercial 
MCMC software package 01- 50,000 samples are obtained for each possible combination of data. Since the resulting 
number density of samples in parameter space is proportional to the joint likelihood, the samples arc then counted 
into bins and those one-dimensional distributions are automatically marginalised over the remaining parameters. 
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TABLE I: Prior for each combination of four data sets: The Las Campanas (L) and IRAS (I) redshift surveys, galaxy peculiar 
velocities (P) and distant supernovae (S) and associated priors. "CDM" means only cold dark matter models were used in the 
likelihood function calculation. "Flat" means that only models of flat universes were considered and the likelihood function 
was generated under that assumption. 



IV. RESULTS 

Table ITT1 lists the best fit values for h, f2 m and as derived from the MCMC samples. When available, the la limits 
are given. "— " means the value was not constrained within the priors. The preferred values for the three parameters 
when all data sets are combined are h = 0.57^®'^, Q m = 0.28 ± 0.05 and a s = 0.87^q'q 5 . This is in the case of a flat 
universe so = 0.72 ±0.05. The best fit values when allowing for non-flat geometries of the universe are given by the 
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combination Las Campanas, peculiar velocities and supernovae; h = 0.56l : 16 , f2 m = 0.3l8; 06 an< ^ a & = 0-87 ± 0.07. 
For that combination, Q A = 0.78±g;?| and fi k = — 0.14±^;^. For convenience, both Q\ and f2k are marginalised over 
in the following pages. Figs. ^ and [3] show the marginalised one-dimensional posterior distribution of the Hubble 
parameter h, the matter density f2 m and the mass density fluctuation parameter as calculated by the MCMC routine. 

Peculiar velocity data seem to prefer higher values of as (panel 3C of figure ^) such that when combined with 
supernovae data, a$ > 1 at the la confidence level (panel 3C of figure 

The combination of the peculiar velocities and supernovae data provides h > 0.6 at la. The combination of 
IRAS with peculiar velocities also prefers 0.55 < h whereas all other combinations prefer lower values of the Hubble 
parameter. The constraints on h obtained here are wide compared to other estimates, e.g. |(|. Table 3 of [l(| shows 
that the latest determinations of h are all in the range preferred here: 0.4 < h < 0.77 at la in the full joint analysis. 

In the combination of Las Campanas and peculiar velocities likelihood functions (panels IB, 3B and 6B of Figure 
respectively), £l m becomes less tightly constrained compared to the individual data sets and shifts to higher values. 
Note that the la intervals for f2 m all overlap and the resulting interval from the combination of all data sets f2 m = 
0.28 ± 0.05 is almost identical with the latest results from [ly , 17 m = 0.29 ± 0.07. This demonstrates that owing to 
the non-linear relation between cosmological parameters, the resulting probability distribution function of derived 
parameters is not simply the product of the distributions obtained individually from each data set. 

Two trends appear in these results that seem related to whether the Las Campanas data are used. The combinations 
that include the Las Campanas data tend to favour lower values of the Hubble parameter and of the mass density 
fluctuation parameter. On the other hand, when IRAS or peculiar velocities data are used higher values of the 
Hubble parameter are preferred, more in accordance with the HST key project or WMAP [l^ results. This is at 
the expense of a higher value of as as shown in panel 6C of Figure [3 



data sets 


h 




0"8 


L 




< 0.11 


Q1+ - 07 
u - Ji -0.08 


I 




o n+ 011 

u - ±o -0.08 


> 0.85 


P 




0.20j™ 


1.48_ .44 


S 




o.28i° : } 5 




L + I 




0.168:8 


n 07+0. 05 

U -°'-0.06 


L + P 


0.4818:2b 


0.3818S 


n Q7+0.06 
U -°'-0.09 


L + S 


> 0.34 


< 0.34 


n Q7+0.08 


I + P 


> 0.58 


n O+0.08 
u - z -0.07 


1.46_ . 59 


I + S 


56+ ' 4 

u - oo -0.19 


0.25 ±0.06 


0.91.Q.05 


P + s 


> 0.63 


0.2818:?! 


1.48.Q.3 


L + I + P 


> 0.33 


o.36i8:i2 


o.85l8:8^ 


L + I + S 


63+ - 21 

U.OO_ 23 


0.26 ±0.05 


0.8618:8! 


L + P + S 


5fi+°- 21 


n q+0.09 
u -°-0.06 


0.87 ±0.07 


I + P + s 


> 0.58 


25 +0 ' 06 

u - ZJ -0.05 


1 0+0.39 

± -°-0.26 


L + I + P + S 


^7+ ' 16 


0.28 ±0.05 


o.87i8:8 5 



TABLE II: Best-fit values for three cosmological parameters from different combinations of four data sets: Las Campanas (L), 
IRAS (I), Peculiar Velocities (P) and Supernovae (S). 



V. CONCLUSIONS 

Four different cosmological sets of data provide constraints on ft, = 0.571° '14, ^ m = 0.28 ± 0.05 and as = 0.871qq5 
under the assumption of a flat universe, and h = 0.56lo^g, fl m = O.3lo;5jg and as = 0.87 ± 0.07 for all geometries. 
The full combination of these data sets provides remarkably tight constraints on f2 m , h and as . The resulting values 
of /Omega m and as are in excellent agreement with parameter estimation presented in |l6(: f2 m = 0.27 ± 0.04, 
as = 0.84 ± 0.04 although the constraint on the Hubble parameter, h = Q.71I5 84, differs somewhat. 
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In general, the three parameters display a different behaviour depending on the data and priors. The Hubble 
parameter is ill-determined in each case and the tightest contraint is consistent at lcr with most h estimates (see |(| 
and references therein). The matter density Q m is fairly well determined in each analysis and all results are in good 
agreement independently of which data are included or which priors are assumed. The amplitude of mass fluctutations 
parameter erg is more difficult to constrain. As in the litterature (see e.g. table 4 in [lfj), several values have been 
derived that are not necessarily consistent and depend strongly on which data are used. Unlike the Hubble parameter 
which is allowed a wide range of values whatever the analysis, observations lead to rather tight constraints on a% but 
different ones, thereby suggesting that the discrepancy is due to errors or misestimates in the data rather than to 
priors. 
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FIG. 1: Constraints on the Hubble parameter h (first column), the matter density fi m (second column) and the mass density 
fluctuation parameter as (third column) obtained from different combinations of large-scale structure data and supernovae. 
Each row corresponds to a different combination. The first four rows are each data set separately. The supernovae likelihood 
function does not constrain <rg. 




FIG. 2: Constraints on the Hubble parameter h (first column), the matter density fi m (second column) and the mass density 
fluctuation parameter as (third column) obtained from different combinations of large-scale structure data and supernovae. 



